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Abstract.  Unsteady  expansion  of  the  vapour  into  vacuum  from  the  flat  surface  was  studied  by  the  direct  simulation 
Monte-Carlo  method.  The  vapour  was  considered  as  monatomic  gas  (the  model  of  hard  spheres).  The  flow  stmcture  and 
parameters  in  the  ranges  of  flow  regimes  from  free-molecular  to  continual  were  investigated. 


INTRODUCTION 

Unsteady  outflow  of  the  vapour  (gas)  into  vacuum  from  the  flat  surface  is  related  to  a  number  of  classic  problems 
of  the  rarefied  gas  dynamics.  For  the  collisionless  gas  expansion  this  problem  has  the  self-similar  analytical  solution 
[1].  The  plane  unsteady  expansion  of  the  gas  into  vacuum  (the  plane  rarefaction  wave)  is  the  continual  analogous  to 
this  flow.  For  the  indicated  flow  in  case  of  the  perfect  gas  the  self-similar  analytical  solution  is  also  known  [2]. 
Generally  the  expansion  of  the  monatomic  gas  into  vacuum  in  the  presence  of  collisions  is  not  self-similar  and  the 
analytical  solution  in  that  case  is  missing.  The  problem  C2in  be  solved  by  approximate  approaches  or  numerical 
methods.  It  was  carried  out  a  sufficient  amount  of  this  kind  of  investigations.  An  important  contribution  on  solution 
of  this  problem  was  made  by  Cercignani  [3],  Anisimov  [4],  Sibold  and  Urbassek  [5,  6],  Kelly  [7]  and  a  number  of 
other  scientists.  In  the  framework  of  these  researches  a  great  attention  was  given  to  the  flow  features  in  the 
nonequilibrium  Knudsen’s  layer,  the  translational  relaxation  process  and  other  significant  properties  of  the  flow. 
Obtained  by  now  results  are  commonly  used  in  various  applications  related  to  the  unsteady  expansion  of  the  vapour 
into  vacuum  (e.g.  pulsed  laser  ablation  of  materials).  At  the  same  time  the  present  results  of  computational 
investigations  are  not  exhaustive  on  the  range  of  analyzed  conditions.  It  is  also  required  more  detailed  analysis  of 
common  mechanisms  of  a  space-time  evolution  of  a  flow  stmcture.  In  practice  of  application  of  the  plane  unsteady 
exp2insion  of  the  gas  into  vacuum  model  for  the  solution  of  more  complex  real  tasks  with  nonplanar  geometry  the 
mistakes  concerned  with  the  incorrect  understanding  of  a  number  of  properties  of  this  flow  are  often  taken  place.  In 
this  report  an  attempt  to  consider  this  problem  in  detail  and  in  corpore  was  made. 


NUMERICAL  METHOD  AND  PROBLEM  STATEMENT 


For  the  simulation  the  DSMC  method  is  used  (Berd’s  version  [1],  NTC  —  model  of  hard  spheres).  The 
computational  domain  includes  a  space  between  an  evaporating  surface  and  an  external  boundary  placed  at  a 
sufficient  distance  from  the  evaporating  surface  L.  Initial  and  boundary  conditions  are  formulated  as  follows.  At  the 
initial  time  moment  (t  =  0)  there  are  no  particles  in  the  computational  domain  (x  >  0).  At  t  =  0  the  gas  begins  to 
outflow  from  the  evaporating  surface  {x  =  0).  The  velocity  distribution  function  of  p2irticles  evaporating  from  the 
surface  at  r  >  0  is  supposed  to  be  constant  (independent  on  time)  and  semi-maxwellian 
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Here  is  the  equilibrium  vapour  concentration,  7",^.  is  the  surface  temperature,  are  the  velocity 

components  of  particles,  R  is  the  gas  constant.  Particles  turning  back  to  the  source  surface  are  excluded  from  the 
simulation  process  that  corresponds  with  the  condition  of  their  total  condensation.  Particles  reached  the  external 
boundary  of  the  computational  domain  x  =  L  are  also  excluded  from  the  simulation  (the  external  boundary  is 
supposed  to  be  the  condensing  surface).  The  modelling  of  the  gas  evaporation  into  vacuum  at  such  conditions  is 
correct  for  a  hypersonic  flow  at  the  external  boundary.  This  requirement  was  always  fulfilled  in  simulations.  For  the 
decreasing  of  the  amount  of  computations  under  the  realization  of  DSMC  method  the  distribution  function  (1)  is 
simulated  in  the  limited  range  of  absolute  values  of  the  velocity.  The  maximum  value  of  velocity  components  under 

investigation  was  =4/^,  P  (2RT .  The  modelling  time  corresponded  to  T  =  tlt^  <10^  and  the 

location  of  the  external  boundary  of  the  computational  domain  corresponded  to  L  <10*  and  —  mean  time 

and  mean  free  path  defined  by  equilibrium  parameters  near  the  surface).  For  the  simulation  the  parallel  DSMC 
algorithm  [8]  was  used. 


MODELLING  RESULTS  AND  THEIR  ANALYSIS 


Initial  Stage  Of  Outflow.  Free-Molecular  And  Nearly  Free-Molecular  Flow  Stages 

In  Fig.  1  the  self-similar  profiles  (dependences  on  ^  =  fi>c / 1 )  of  relative  parameters:  velocity  /h  ,  density 
.  longitudinal  T  ^  IT^,  transverse  T,.  /T,,.  and  total  T  IT^.  temperatures,  flux  of  particles  0un/ and  Mach 

number  M  for  the  collisionless  flow  as  well  as  the  results  of  calculation  of  these  parameters  at  r  =  10”^,  10”'  and  1 
are  given.  Results  obtained  practically  coincide  with  each  other  and  with  self-similarly  profiles  for  the  collisionless 
gas  expansion.  Reliable  design  data  are  limited  by  following  range  <  3.2  —  3.4  .  Its  right  boundary  is  marked 
with  figure  1  (Fig.  1).  The  flow  on  this  stage  is  characterized  by  the  rapid  density  decreasing  and  the  nearly  linear 
increasing  of  gas  velocity  while  moving  away  from  the  surface.  The  dashed  line  2  in  Fig.  1  corresponds  to  asymptotic 
profiles  of  relative  parameters  at  great  values  of  ^ .  The  flow  is  completely  translationally  nonequilibrium.  The 
transverse  temperature  T,,  is  constant.  The  longitudinal  temperature  T ^  is  monotonously  decreasing  while  moving 
away  from  the  surface.  The  sonic  point  (M  =  1)  corresponds  to  ^  =  0,337.  Parameters  at  the  sonic  point  are  the 
following:  tki  =  0,795  ,  n/n^  =0,317  ,  T  /r„,  =0,757,  fiunln^=  0,252  . 

At  r  =  1  the  averaged  number  of  collisions  of  particles  under  the  outflow  into  vacuum  is  about  0. 1 .  The  value 
r  =  1  may  be  assumed  as  the  boundary  of  the  free-molecular  stage  of  the  gas  outflow.  At  r  =  3  the  number  of 
collisions  is  increasing  up  to  1 .  The  influence  of  collisions  on  profiles  of  gasdynamic  parameters  becomes  more 
appreciable  but  the  difference  of  design  profiles  from  the  self-similar  profiles  remains  insufficient  (less  5  -  1 0  %  for 
T ^  and  T ^  ,  less  2  —  3  %  for  «  and  n  ).  The  rtinge  1  <  r  <  3  may  be  assumed  as  the  stage  of  the  nearly  free- 

molecular  regime.  At  this  stage  the  transition  from  collisionless  to  collisional  flow  regimes  occurs. 


a  b 

FIGURE  1.  Self-similar  profiles  fiu,  n/n^,  pun/n^  M  (a)  and  T„  T„  T  (b)  at  collisionless  flow  and  design  results  at  r=10'^  (•),  10'' 
(□),  1  (o);  1  —  the  boundary  of  the  reliable  design  data  region,  2  —  asymptotic  profiles  of  parameters  at  great  values  of  ^ 
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Transitional  And  Continual  Stages  Of  Flow  Regimes 


Let's  consider  the  flow  evolution  at  r  >  3  .  In  Fig.  2  profiles  pu  ,  n! and  pun  / on  coordinate  x  ! (left) 
and  self-similar  coordinate  ^  (right)  at  r  =  3—10  000  are  given.  There  is  also  presented  the  self-similar  solution  for 
the  collisionless  flow  (curve  0)  in  Fig.  2  (right).  The  origin  of  design  curves  in  Fig.  2  corresponds  to  the  middle  of 
the  first  cell  (  12  ).  With  the  T  growth  the  flow  suffers  fundamental  changes.  It  becomes  non-self-similar  as  a 

whole.  This  is  connected  with  the  appearance  of  the  characteristic  scale  (mean  free  path)  defining  sizes  of  zones  of 
the  translational  relaxation.  Profiles  pu  and  are  not  self-similar  near  the  surface  but  they  are  close  to  self¬ 

similar  profiles  at  sufficiently  great  Sizes  of  the  region  of  the  approximate  self-similarity  of  pu  and  nln^  are 
growing  with  the  T  increase.  At  r>1000  the  approximate  self-similarity  of  these  parameters  takes  place  at 
^  >0.1.  Profiles  pun  /  n^,  are  approximately  self-similar  for  the  whole  flow  region. 

With  the  r  increase  the  number  of  collisions  of  particles  is  growing  in  the  gas,  the  reverse  flow  of  particles  to 
the  surface  appearing  and  growing.  This  results  in  decreasing  pu  and  increasing  near  the  surface.  In  some 

zone  near  the  surface  the  magnitude  pun! n^  is  changed  insufficiently  (Fig.  2e,  f).  At  r  >  100  the  value  of  pun! 
near  the  surface  remains  approximately  constant  and  equal  to  0.24.  The  10  %  decrease  in  punt n^  occurs  at  the  distance 
from  the  evaporating  surface  x  !  A^=  (0,3  —  0,35)r  or  ^  =  0,33  -  0,4  . 


f 


FIGURE  2.  Dependences  Pu  ,  Yl  I  and  pin  I  on  X  /  {a,  c,  e)  and  ^  {h,  d,f)  at  various  time: 

r  =  0  (curve  0),  3  (I),  10  (2),  30  (i),  100  {4),  300  (5),  1000  ((5),  3000  (7),  10000  (,S) 
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The  governing  regularities  of  this  flow  are  connected  with  the  process  of  the  translational  relaxation.  In  Fig.  3 
profiles  IT and  T ^  IT at  r  =  30  {a)  and  300  (h)  are  given.  They  represent  the  course  of  the  translational 

relaxation  with  the  r  growth.  The  horizontal  line  in  Fig.  3  corresponds  to  the  value  of  T,.  /r„.  at  x  — >  oo  .  The 
current  position  of  the  gas  leading  edge  is  denoted  as  =4f/^.  At  r  =  30  the  traslational  equilibrium  isn’t  settled 
due  to  the  insufficient  number  of  collisions.  At  r  >  100  the  zone  of  the  quasi-equilibrium  flow  is  formed  at  some 
distance  from  the  surface.  Within  the  bounds  of  this  zone  values  of  T^.  and  T are  close,  <  7',. . 

The  following  condition  can  be  assumed  as  a  criterion  of  the  traslational  equilibrium 


Within  the  bounds  of  the  Knudsen's  layer  S  >  S  ,  a\.  it’s  external  boundary  5  =  S  .\xi  Fig.  4  data  on  the  position  of 
the  external  boundary  of  the  Knudsen’s  layer  at  different  r  for  e  =  0,05  are  given.  At  first  at  r  >  100  the 

Knudsen’s  layer  thickness  sufficiently  rapidly  decreases  and  then  at  r  >  1000  it  is  changed  slightly  and 

approximately  equal  to  (II  —  As  the  time  passes  parameters  at  the  point  Xj^^  are  changed:  values  of  nln^ 

and  T  IT ^  increase,  pu  and  M  decrease.  M  number  at  the  external  boundary  of  the  Knudsen’s  layer  in  the  range  of 
r  =  10^  -  lO"*  is  changed  from  1.48  to  0.775. 

The  governing  regularities  of  the  translational  relaxation  beyond  the  bounds  of  the  Knudsen’s  layer  are  the 
following.  The  magnitude  T^  monotonously  decreases  and  the  curve  7"^,  has  its  minimum.  The  minimum's 

coordinate  x^  is  moving  away  from  the  surface  with  the  r  growth.  The  dependence  x^  (r)  (Fig.  4)  at  r  >  30  can 
be  approximated  by 

X  IX  =0,67r‘-''',  B  =  1,12^  =  0,75r‘’•’^  (3) 

me’  ’  “m  ,-3  ''  ' 


FIGURE  3.  Dependences  /  Tfy  (a,  b>,  /  T^^■  (a,  b),  T  /  T^^,  (c)  and  M  (d)  on  X  !  at  various  time:  a  -  r=30,  b  -  r=300, 

candd-  i=3  {!),  10(2),  30  (3),  100  (4),  300  (5),  1000  (d),  3000  (7),  10000  (S) 

At  X  >  x^  the  value  of  T rapidly  increases  approaching  to  the  value  of  IT^  =  1  corresponding  with  the  free- 

molecular  mode  of  the  flow.  Such  behavior  of  curves  T is  accounted  for  particles  practically  collisionless  moving 
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in  the  region  adjoining  to  the  gas  leading  edge.  This  feature  of  changing  T^,  at  x>x^  can  be  used  for  the 
determination  of  the  boundary  of  the  flow  transition  to  the  free-molecular  mode  of  the  expansion.  As  such  boundary 
we  shall  consider  die  intersection  of  the  design  curve  T^.  linearly  extended  (the  dashed  line  in  Fig.  3a)  with  the  line 

T^,  /T^,  =  1 .  Let's  denote  this  point  as  x ^  .  The  dependence  x  (r)  is  given  in  Fig.  4.  Design  data  at  r  >  30  are 
approximated  by; 

=  (2.9  -  3^)r,  3,5  .  (4) 

Regularities  of  and  T^,  changes  identically  determine  profiles  T.  The  typical  feature  of  the  T  profiles  is  the 
presence  of  the  minimum  at  the  point  and  the  tendency  of  T  to  the  value  T  =  0,677’^^,  in  the  region  of  the  free- 
molecular  expansion  (X  >  X  ).  Profiles  u  and  T  define  Mach  number  profiles.  The  presence  of  the  T  minimum 
results  in  appearing  the  local  maximum  of  the  M  profile  at  the  point  x^  at  r  >  1 00  (Fig.  3tf).  M  number  at  the  local 
maximum  point  is  increasing  with  the  r  growth.  Simulations  carried  out  describe  the  flow  in  the  region 
^<3,2  — 3,4.  It  was  shown  above  M~x/t  for  the  free-molecular  region  of  the  gas  expansion.  The  hypothetical 
form  of  the  M  profile  at  x  >  for  r  =  3000  is  given  in  Fig.  3ct  (dotted  line  9). 

At  r>100  the  form  of  gasdynamic  parameters  profiles  (Fig.  2  and  3)  points  out  the  formation  of  two  flow 
regions  with  the  distinction  in  kind  of  regularities  of  parameters  changes:  the  subsonic  quasi-steady  region  with  the 
thickness  A.  and  the  supersonic  flow  region  downstream  located  (Fig.  4).  Both  regions  have  the  fairly  complex 
internal  structure.  In  Fig.  4  data  on  the  sonic  point  location  (M  =  1)  are  presented  while  in  Fig.  8  parameters  at  this 
point  are  given.  The  coordinate  of  the  sonic  point  .x.  is  monotonously  increasing  with  the  r  growth.  At  r  >  30 
dependences  x,/A^  and  on  r  can  be  approximated  by 

x./;i,  (5) 

At  r  =  =  200  (Fig.  4)  values  of  x,^  and  x,  are  equal.  At  r  >  the  external  boundary  of  the  Knudsen’s  layer 

is  located  in  the  region  of  the  subsonic  flow.  In  the  range  100<r<r^  wehavex^j^^  >  x. .  Within  the  bounds  of  this 
short  interval  r  the  flow  at  the  sonic  point  is  nonequilibrium.  At  the  sonic  point  at  r  =  30 -10  000  (Fig.  5) 
parameters  n,  /  n^,  pu,  and  T-  IT^y  are  changed  slightly  (especially  at  r  >  1000  ).  This  feature  can  be  used  for  the 
approximate  statement  of  boundary  conditions  at  the  surface  under  the  present  problem  solution  by  the  continual 
gasdynamic  approaches  (Euler  or  Navier  —  Stokes  equations). 


FIGURE  4.  Qa.sdynamic  structure  of  the  flow  at  transitional  and 

...  ^  _  FIGURE  5.  Parameters  at  the  sonic  point  of  the  flow 

continual  stages  at  coordinates  X  /  ,  T 

The  flow  in  the  subsonic  region  is  characterized  by  the  fairly  rapid  decrease  of  gradients  of  gasdynamic 
parameters  with  the  x  increase  (Fig.  2  and  3).  With  the  r  growth  this  tendency  is  gaining  strength.  At  the  stage  of 
the  developed  continual  regime  (  r  >  1000  )  in  the  subsonic  region  the  zone  of  the  rapid  acceleration  of  the  gas  up  to 
M  =  0,8  -  0,9  with  die  extension  (20  —  30)  T,  and  the  transonic  zone  with  the  much  great  extension  can  be  pointed 
out.  The  size  of  the  transonic  zone  is  growing  with  the  T  increase.  Downstream  the  supersonic  unsteady  region  is 
located.  Within  the  bounds  of  this  region  the  rapid  growth  of  /kt  ,  decrease  of  n!  and  pun  /  are  observed. 

At  r  >  100  downstream  from  the  Knudsen's  layer  the  region  of  the  continual  (quasi-equilibrium)  flow  is  located. 
The  extension  of  this  zone  is  growing  with  the  T  increase.  As  an  external  boundary  of  the  continual  region  can  be 
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assumed  the  coordinate  X^.  where  condition  (2)  is  true.  The  location  at  £  =  0,05  is  presented  in  Fig.  4.  The 

extension  of  the  region  of  the  continual  flow  is  equal  to  .  At  some  magnitude  r  =  the  external 

boundary  of  the  region  of  the  continual  flow  is  confluent  with  the  external  boundary  of  the  Knudsen's  layer 
(Fig.  4).  At  £■  =  0,05  we  have  =  100  .At  T  <T^  the  Knudsen’s  layer  as  a  zone  of  the  relaxation  of  the  initial 
nonequilibrium  cannot  be  assumed  since  the  flow  located  left  from  the  point  «  is  nonequilibrium  in  the  whole  region. 
Above  the  line  x  the  region  of  the  free-molecular  flow  is  located.  Between  the  lines  x  and  x^  the  region  of 

the  transitional  regime  of  the  flow  is  located.  Its  size  =  x  .  —  x ^  is  growing  with  the  T  increase. 


GENERALIZED  GASDYNAMIC  STRUCTURE  AT  COORDINATES  r 


The  summarizing  of  design  data  allows  to  represent  the  time  evolution  of  the  gasdynamic  structure  for  the  flow 
as  a  whole  (Fig.  6).  At  coordinates  r  the  whole  flow  field  can  be  divided  into  three  regions;  free-molecular  F, 
transitional  T  and  continual  C  flows.  At  r  >  3  the  region  of  the  free-molecular  regime  has  two  subregions:  the  first 
one  is  located  above  the  line  rf  corresponding  to  the  boundary  of  the  free-molecular  regime  ^ ,  the  second  one  is 


placed  below  the  line  rp  corresponding  to  the  boundary  of  the  vapour  layer  with  the  thickness  X  = 

adjoining  to  the  evaporating  surface.  Left  from  the  point  r  (t  <3)  the  flow  is  completely  free-molecular  or  nearly 
free-molecular.  The  region  of  the  transitional  regime  is  located  between  the  lines  prf  and  knc.  The  line  knc  divides 
the  transitional  flow  from  the  continual  one.  It  consists  of  the  external  boundary  of  the  Knudsen’s  layer  and  the 


distant  boundary  of  the  region  of  the  continual  flow  .  Right  from  the  line  knc  the  region  C  is  located. 


The  time  evolution  of  the  gasdynamic  structure  of 
the  flow  has  four  typical  stages  (Fig.  6):  1  —  free- 
molecular  ( r  <  1 ),  2  —  nearly  free-molecular 

(l<r<3),  i  —  transitional  (3<r<100)  and  4  — 
continual  (r>10()).  Here  and  above  names  of  flow 
stages  correspond  to  names  of  mostly  dense  flow 
regimes  realized  at  the  corresponding  stage.  At  stages  1 
and  2  the  flow  regime  F  takes  place.  At  stage  i  the  flow 
includes  the  region  of  the  transitional  regime  T 
<^f„)and  two  zones  of  the  region  F  (^  <  4^^, 


FIGURE  6.  Generalized  gasdynamic  structure  of  the  flow  at 
coordinates  4  >  ^ 


and  4  >  )•  At  the  forth  (continual)  stage  the  flow 

includes  regions  with  three  flow  regimes. 


The  Knudsen’s  layer  (  ^  <  4f^n  )  including  the  free-molecular  sublayer  ( ^  <  4 the  zone  of  the  transitional 
flow  (4fK  <  ^  ^Kn  ^  placed  near  the  surface.  Farther  the  regions  of  continual  <4<^^),  transitional 


(4^  <  4  )  3nd  free-molecular  (4  >  )  flows  are  located. 
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